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TRODUCTION 


A  propeller  operating  in  the  turbulent  wake  of  an  axisymmetric  body  with  appendages 
encounters  wake  non- uniformities  which  result  in  spatial  and  temporal  fluctuations  of  blade  angle-of- 
attack.  These  angle-of-attack  fluctuations  result  in  unsteady  blade  loadings  and  the  generation  of 
propeller  noise,  and  the  noise  sources  are  characterized  by  three  types  of  unsteady  force  mechanisms: 
a)  turbulence  injestion;  b)  vortex  shedding;  and  c)  blade-rate.  The  first  two  mechanisms  typically 
generate  continuous  spectrum  (broadband)  radiated  noise  levels  while  blade-rate  forces  generate 
discrete  frequency  noise  levels  at  various  blade-passage  frequencies  and  harmonics.  This  paper  will 
address  the  reduction  of  blade-rate  noise.  [A  .ti  5  The  hull  boundary  layer  behind  an  appendage  (i.e. 
figure  l )  is  characterized  by  a  complex  velocity  field  typically  having  velocity  excesses  at  inner  radii 
and  velocity  defects  at  the  outer  radii.  This  type  of  velocity  field  has  a  complex  harmonic  content 
distribution  and  its  effect  on  blade-rate  noire  cannot  be  predicted  without  detailed  examination  of 
the  wake  and  the  radial  distribution  of  propeller  blade  forces.  The  reduction  of  blade-  rate  noise  and 
thus  vehicle  vibration  provides  the  motivation  for  the  application  of  skew  in  propeller  design. 


Techniques  are  available  to.  computing  unsteady  forces  and  skew  distributions,  and  these 
methods  range  from  low-aspect  ratio  approximations  to  unsteady  airfoil  theory  to  complete  unsteady, 
lifting-surface  methods.  However,  since,  no  method  was  available  to  systematically  determine  an 
optimum  skew  distribution,  the  propeller  skew  optimization  program  SKEWOPT  (Gr^enblatt,  1978;  ^ 
and  Parsons  and  Greenblatt,  1978)  was  developed.  SKEWOPT  determines  a  quadratic  or  cubic  skew 
distribution  using  an  optimization  technique  which  finds  the  set  of  parameters  for  which  a  user- 
defined  linear  combination  of  the  unsteady  force  and  moment  amplitudes  are  minimized.  SKEWOPT 
was  written  for  use  in  ship  propeller  design;  however,  the  intention  was  that,  after  minor  modifica- 
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tions  such  as  the  inclusion  of  higher  order  harmonic  groups,  it  would  be  suitable  for  torpedo  and 
submarine  propeller  design.  Since  the  force  calculation  method  in  SKEWOPT  was  not  sufficiently 
documented,  the  method  was  replaced,  and,  due  to  the  difficulties  encountered  in  modifying 
SKEWOFT,  a  new  progmm,  patterned  after  SKEWOPT,  was  written. 

THE  SKEW  OPTIMIZATION  PROBLEM 

VELOCITY  FIELD 


Figure  1  shows  one  quadrant  of  a  wake  behind  an  axisymmetric  body  having  four  identi¬ 
cal  control  surfaces.  The  spatial  variations  in  the  circumferential  and  radial  directions  are  due  to  both 
potential  and  viscous  effects,  and  it  is  this  type  of  velocity  field  that  plays  an  important  part  in  the 
design  of  wake-adapted  propellers.  While  circumferentially  averaged  velocity  profiles  can  be  used  for 
propeller  design  calculations,  the  determination  of  a  skew  distribution,  for  the  reduction  of  unsteady 
forces,  requires  that  the  spatial  variations  of  the  inflow  velocity  field  be  considered.  Since  the  spatial 
velocity  distributions  are  periodic  and  continuous,  they  may  be  represented  in  terms  of  a  Fourier 
series  or  a  complex  exponential.  For  example,  the  axial  component  of  the  velocity  at  a  position  ( r,B ) 
can  be  expressed  as 


js&a. 

U„ 


«o(r) 

2 


JJ  [a*(r)cos(n0)  +  b„(r)  sin(nJ)  |  -  Re 


a.  oo 


c„  Cto* 


(1) 


where  U„  is  the  free  stream  velocity,  a„(r) ,  a^r) ,  b*(r)  and  c*(r)«Mr)'-ib*(r)  are  the  Fourier  coeffi¬ 
cients  and  Re  ( )  denotes  the  real  part  (only  the  real  part  will  be  used).  Since,  the  term  a*  does  not 
vary  in  6 ,  it  is  associated  with  the  steady  state  thrust  and  torque,  and  the  additional  terms  are 
sinusoidal  fluctuations  of  the  inflow  velocity  which  produce  the  unsteady  forces  and  moments. 


Once  thr.  velocity  field  and  propeller  geometry  arc  determined,  the  problem  of  determin¬ 
ing  an  optimum  skew  distribution  requires  the  formulation  of  a  nonlinear  programming  problem 
which  includes  an  unsteady  force  calculation  method.  First,  the  calculation  of  the  unsteady  forces 
will  be  considered. 
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UNSTEADY  FORCE  CALCULATION  METHOD 


The  original  version  of  SKEWOPT  (Greenblatt,  1978)  had  both  a  two-dimensional, 
unsteady  and  a  more  time  consuming  lifting-line  method  available  to  calculate  blade  forces.  To  over¬ 
come  computational  problems,  the  original  SKEWOPT  force  calculation  methods  were  replaced  by  a 
method  developed  by  Thompson  (1976).  Briefly,  his  method  divides  the  propeller  blade  into  strips 
which  are  considered  two-dimensional  airfoils.  Included  in  the  method  are  a)  the  two-dimensional 
unsteady  airfoil  theories  of  Sears  (1941)  and  Horlock  (1968)  which  allow  consideration  of  sinusoidal 
velocity  fluctuations  normal  and  parallel  to  the  inflow  velocity,  and  b)  corrections  to  the  blade  lilt 
force  due  to  the  presence  of  adjacent  propeller  blades.  The  effect  of  blade  camber  (from  Naumann 
and  Yeh,  1973)  has  also  been  added. 

In  order  to  reduce  blade-rate  propeller  forces,  six  unsteady  forces  and  moments  must  be 
considered.  As  shown  in  figure  2,  they  are  the  thrust  FB,  the  side  forces  Fx  and  F,,  the  torque  Tt,  and 
the  bending  moments  Tx  and  Tr  For  an  axisymmetric  body,  one  is  not  concerned  with  differentiat¬ 
ing  between  the  x  and  y  components  of  the  side  forces  and  bending  moments;  therefore,  the  two  side 
forces  and  bending  moments  were  combined  to  yield  the  maximum  side  force,  F„  and  maximum 
bending  moment,  T,.  The  derivation  of  the  force  calculation  method  and  extensions  are  given  in 
Thompson  (1976)  ana  Mautner  and  Blaisdell  (1987). 

It  has  been  shown  (Thompson,  1976)  that  only  certain  inflow  harmonics  contribute  to  the 
unsteady  forces  and  moments.  They  are  the  harmonics  for  which  the  order  is  some  multiple  (m)  of 
the  number  of  blades  (Nb).  Harmonics  of  order  mNb,  gives  rise  to  an  unsteady  force  and  moment  in 
the  z  direction,  Fa  and  M,  ,  while  the  side  forces,  Fs  and  Fr  and  bending  moments,  M*  and  My,  are 
generated  by  inflow  harmonics  of  order  mNb«±l.  The  equations  for  the  thrust,  torque  and  maximum 
side  force  and  bending  moment,  for  a  given  harmonic  group  m,  are 
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(2) 


F.W“Nb  g  [  L,{lBNk) 

T.W-Nh  £  f  1  r,  tan(ft)  c*"*  I  (3) 

j>k  l  J*  I 

♦ 

F.<1-&nJ  £  [  ]  _  WK4)  e***  +  ttN,  £  [  L,,.*,.,,  ]  _  e-"*  (4) 

T.'-'-HiN*  £  [  ],  r,  *  UN  £  [  I,*,,.,,  ]_  r,  (5) 

where  L  is  the  lift  force,  ^  is  the  skew,  P  is  the  number  of  blade  strips,  fi  is  the  pitch  angle  and  the 
index  m  includes  harmonic  groups  m- 1,2, 3,4.  It  should  be  noted  that  the  contributions  of  the  dif¬ 
ferent  harmonic  groups  are  given  separately  since  the  unsteady  forces  due  to  different  harmonic 
groups  fluctuate  at  the  frequencies  mN^O  where  ft  is  the  propeller  rotation  rate. 

SKEW  OPTIMIZATION  MODEL 

To  determine  the  optimum  skew  distribution  die  above  force  calculation  method  has  been 
incorporated  into  a  nonlinear  programming  problem.  Due  to  die  fact  that,  in  general,  all  forces  cannot 
be  minimized  simultaneously,  a  scalar  cost  function  formed  from  the  weighted,  linear  combination  of 
the  forces  and  moments  is  minimized  The  cost  function  Fc  is 


Wj^F.W  Wj^F.W  Wj-frW  wj-frrj-* 
.05F.  +  .05?.  +  .05T,  +  .05T, 


(6) 


where  the  weights,  Wj,  are  normalized  such  that  their  sum  over  both  m  and  j  (■!,. .,4)  equals  1. 


Since  the  Fc  depends  upon  the  skew  distribution,  ^(r),  one  could  solve  for  the  optimum 
skew  via  a  variational  or  n  dimensional  parameter  optimization  technique.  However,  a  more  feasible 
approach,  and  the  one  used  here,  is  to  use  a  few  parameters  in  describing  the  skew  and  perform  an 
optimization  search  in  a  limited  parameter  space.  To  accomplish  this,  the  skew  distribution  is 
represented  by  either  a  cubic  or  quadratic  distribution  having  a  straight  line  section  with  ^-0.  The 


distributions  (see  Geenblttt,  1978  and  Mautner  and  Blaisdell,  1987)  are 


^tr)-aiJ*br*+cr-fd  er  br^+cr-fd  r,  £  r  £  r» 

*r)-0  r^r$  r.  (7) 

•  * 

where  is  the  hub  radius,  r»  is  the  tip  radius  and  r,  is  the  radius  at  which  the  polynomial  skew  dis¬ 
tribution  starts.  This  type  of  skew  distribution  is  illvistrated  in  figure  3. 

The  skew  distribution  given  by  equation  (7)  has  five  free  parameters  (a,  b,  c,  d,  r,);  how¬ 
ever,  it  is  more  meaningful  to  the  ptxpeller  designer  to  use  parameters  that  have  physical  meaning 
instead  of  using  these  polynomial  coefficients.  The  parameters  chosen  for  use  in  die  current  method 
are  the  skew  at  the  propeller  tip  dp  the  starting  skew  slope  S,  -  dlO*  the  starting  radius  r„  and  die 
skew  slope  at  the  tip  S*  Additionally,  the  physical  restriction  that  ^(r,)*0  is  made  so  that  the  skew 
distribution  is  continuous  thereby  reducing  die  number  of  free  parameters  by  one.  Another  neces¬ 
sary  restriction  is  that  the  skew  distribution  be  smooth  which  implies  that  S,  -  d^rJ/dr-O  if  r,  >  rh. 
If  there  is  no  straight  line  section  r,  -  rk,  Sp  is  not  restricted  but  r.  is  fixed.  In  either  case,  the 
number  of  free  parameters  is  reduced  to  three  for  a  cubic  distribution  and  two  for  a  quadratic  distri¬ 
bution  and  results  in  the  four  possible  skew  distribution  models  presented  in  table  1.  Specification  of 
these  models  allows  the  optimization  search  to  be  carried  out  in  either  a  two  or  three  dimensional 
space. 

The  propeller  skew  distribution  problem  has  now  been  formulated  as  an  optimization 
problem  in  terms  of  a  few  geometric  parameters.  In  order  to  obtain  a  feasible  propeller  geometry,  it 
is  necessary  to  place  some  restrictions,  such  as  a  maximum  allowable  tip  skew,  on  the  geometric 
parameters.  In  doing  so,  the  design  problem  becomes  a  constrained,  nonlinear  optimization  problem 
where  the  optimization  search  is  restricted  to  finding  the  set  of  parameters  which  minimizes  the  cost 
function  Fe  while  satisfying  all  of  the  constraints  placed  on  the  problem.  The  constraints  used  arc 
given  in  table  2  and  are  checked  for  violation  at  a  given  point  in  the  parameter  space.  The  additional 
constraint  that  the  propeller  blade  should  rot  curve  forward  has  been  incorporated. 
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The  constrained,  nonlinear  optimization  problem  is  represented  by 

min  F«(2&)  subject  to  Oj  Q£)  £  0  (8) 

where  JCita  vector  in  parameter  space  which  determines  the  skew  distribution  and  Gj  (X)  describes 
the  constraints.  There  are  many  techniques  for  solving  the  unconstrained  minimisation  problem  (see 
Parsons,  197S);  however,  only  a  few  methods  attack  the  constrained  problem  directly.  One  useful 
technique  is  to  convert  the  constrained  problem  into  an  unconstrained  problem  and  then  use  an 
unconstrained  optimization  method.  This  can  be  accomplished  with  the  use  of  an  external  penalty 
function  which  is  added  to  a  cost  function  whenever  a  constraint  is  violated  (Le.  Gj(X)  <  0  ).  The 
penalty  function  can  be  expressed  as 

P(XA)  -  F(X)  -  r*  £  min  [Gj(X),  0]  (9) 

l 

and  then  an  unconstrained  optimization  technique  can  then  be  applied  to  P (X  ,rk).  If  no  constraint  is 
violated  no  penalty  is  added  and  the  penalty  (Unction  is  the  same  as  the  cost  (Unction.  Since  the 
penalty  added  is  proportional  to  the  constraint  violation,  the  optimization  method  should  be  forced 
towards  a  feasible  region  where  no  constraints  are  violated.  This  will  be  the  case  as  long  as  the  mul¬ 
tiplicative  factor  rk  is  lar^e  enough  («1024).  If  rk  is  too  small,  die  search  may  tend  toward  an  infeasi¬ 
ble  region.  The  optimization  technique  chosen  for  use  in  this  method  is  the  Nelder-Mead  simplex 
search  method  (Nelder  and  Mead,  1963),  and  a  detailed  description  of  the  method  as  coded  in  the 
current  computer  program  is  given  by  Mautner  and  Blaisdell  (1987). 

The  computer  program  which  solves  the  above  nonlinear  programming  problem  is  an 
interactive  program  intended  for  routine  propeller  design  work.  The  program  has  an  interactive  input 
method  to  accept  propeller  data  and  optimization  parameters,  and  the  program  can  be  restarted,  for 
example,  with  different  constraint  values.  In  addition  to  the  optimization  mode,  a  test  mode  is  avail¬ 
able  for  the  calculation  of  unsteady  forces  for  a  given  skew  distribution. 
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SKEW  DISTRIBUTION  CALCULATIONS 


To  demonstrate  the  results  that  can  be  obtained  from  the  currant  skew  optimization  pro* 
gram,  the  forward  propeller  of  a  counterrotating  propeller  set  was  used  to  calculate  several  optimum 
skew  distributions.  The  propeller  was  designed  using  the  velocity  data  measured,  in  a  wind  tunnel, 
by  Nelson  and  Fogarty  (1977)  and  the  propeller  design  programs  of  Nelson  (1972,  1973).  The  operat¬ 
ing  parameters  and  the  unskewed  blade  geometry  are  given  in  table  3  and  schematically  shown  in  fig¬ 
ure  2. 


It  has  been  mentioned  that  only  certain  harmonics  of  the  inflow  velocity  field  contribute  to 
the  unsteady  forces  and  moments.  Since  the  current  propeller  has  six  blades  and  is  operating  in  a 
four  cycle  wake,  only  those  harmonies  which  are  integer  multiples  of  the  blade  number  H,  need  be 
considered,  and  the  forces  of  interest  are  the  unsteady  thrust,  F„  and  torque,  Ta.  To  provide  the  final 
set  of  data  required  for  the  skew  calculations,  a  Fourier  analysis  of  the  input  wake  (figure  i)  was 
performed.  The  results  show,  for  the  six  bladed  propeller,  the  dominance  of  the  12th  harmonic  and 
the  rapid  approach  to  a  nearly  zero  magnitude  of  the  24th  and  higher  harmonics  (Mautner,  1987). 
The  radial  distribution  of  the  12th  harmonic  is  plotted  in  figure  4. 

The  propeller  geometry  and  the  radial  distribution  of  the  12th  and  24th  harmonics  were 
used  in  determining  the  skew  distributions  given  in  figure  5  and  the  magnitude  of  the  forces  and 
moments  presented  in  table  4.  To  provide  a  reference,  the  magnitude  of  the  forces  and  moments  for 
the  unskewed  propeller  were  calculated.  Next,  each  of  the  four  skew  distribution  models  with  the 
appropriate  constraints  (see  tables  2  and  4)  was  specified.  The  actual  values  of  the  constraints  and 
the  resulting  values  of  the  total  forces  and  moments  are  given  in  table  4,  and  the  calculated  skew  dis¬ 
tributions  are  plotted  in  figure  5. 

It  is  clear  from  an  examination  of  the  magnitudes  of  the  forces  and  moments  given  in  tabic 
4  that  a  significant  reduction  of  the  the  total  thrust  and  torque  on  the  propeller  was  achieved  irre- 
gardlcss  of  the  skew  model  used.  It  can  also  be  seen  that  the  degree  of  force  reduction  and  the  blade 
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tip  skew  varies  with  skew  model  with  the  quadratic  #1  model  product  the  lowest  force  and 
moment  values  (  >tO%  reduction  )  and  the  maximum  tip  skew  (  54*  ).  The  radial  distribution  of  skew 
shown  In  flfure  5  Indicates  dose  similarity  in  the  results  of  each  model  and  skew  distributions  which 
will  yield  satisfactory  propeller  |eomc tries  without  severe  restrictions  being  placed  on  the  skew 
model  constraints.  Clearly,  the  choice  of  which  model  to  use  or  what  level  of  force  reduction,  max* 
imum  tip  skew  or  propeller  geometry  is  acceptable  must  be  determined  by  die  application.  It  should 
be  mentioned  that  the  actual  magnitude  of  the  force  or  moment  is  not  as  important  as  die  relative 
force  reduction  from  die  unskewed  to  skewed  propeller  geometries. 

CONCLUSION 

Based  on  SKEWOPT,  enhanced  skew  optimisation  technique  was  developed.  The 
current  method  provides  a  fast  and  efficient  way  to  determine  a  variety  of  cubic  or  quadratic  skew 
distributions  which  minimise  the  unsteady  forces  produced  by  the  various  harmonic  components  of 
the  input  wake.  The  method  has  been  extended  to  include  higher  order  harmonics,  allows  die  investi¬ 
gation  cf  force  magnitudes  due  to  individual  harmonki  and  the  calculation  of  forces  for  specified 

,  t 

I 

skew  distributions.  Since  the  program  is  interactive  and  involves  relatively  short  computation  times,  it 

t  » 
provides  a  valuable  tool  for  propeller  design. 
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Table  1.  Skew  Distribution  Models. 
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Table  2.  Skew  Distribution  Constraints 
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FIGURE  X.  CIRCUMFERENTIAL  VARIATION  OF  THE  MEASURED  INFLOW  VELOCITY  FIELD 


I 


FIGURE  2.  DESCRIPTION  OF  A  TYPICAL  PROPELLER  AND  ITS  GEOMETRY 
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